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ABSTRACT 

For  inhomogeneous iv  broadened  media  such  as  polydiacetylene -4BCMU  films, 
a  phenomenological  model  is  proposed  to  calculate  the  imaginary  part  of  the 
optical  susceptibilitv  as  a  function  of  the  pump  field  intensity.  It  is 
shown  that  the  model  can  qualitatively  account  for  the  anomalous  behavior  cf 
the  absorption  hole  observed  in  recent  experiments.  Possible  cause  of  the 
anomaly  is  discussed. 


INTRODUCTION 

Studies  of  transient  optical  effects  have  become  important  in  the 
development  of  coherent  optical  spectroscopy  and  fast  response  nonlinear 
optical  materials.  Many  experiments  have  been  carried  out  by  the  pump -probe 
technique  in  which  the  strong  pump  laser  induces  the  nonlinear  optical 
response  of  the  sample  to  be  examined  by  the  weak  probe  laser.  Spectral 
hole  burning  or  bleaching  accompanied  by  a  dynamical  Stark  shift  has  been 
observed  in  such  experiments  on  various  materials  such  as  semiconductor 


1  2 

quantum  wells  and  conjugated  polymers.  A  phase-space  filling  model  has 
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been  proposed  which  provides  a  good  explanation  for  the  hole  burning  and 
exciconic  Stark  effect  in  low- dimensional  semiconductors. 

It  is  well-known  that  polymers  exhibit  giant  nonlinear  optical 
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susceptibilities  with  fast  responses.  They  have  attracted  much  attention 
in  recent  years  for  their  potential  applications  in  future  generations  of 
high-speed  signal  processing  devices.  It  is  therefore  of  great  interest  to 
investigate  the  mechanism  responsible  for  the  nonlinear  optical  behavior. 

In  recent  experiments  on  polydiacetylene  (PDA) -4BCMU, ^  the  hole  width 
measured  in  the  absorption  spectrum  is  found  to  change  with  the  pump  laser 
intensity  in  an  unexpected  fashion.  It  decreases  first  and  then  increases 
as  the  intensity  of  the  incident  beam  increases.  On  the  other  hand,  the 
absorption  depth  increases  slowly  at  first  and  rises  sharply  as  the  pump 
intensity  steadily  increases.  Vie  attempt  in  this  paper  to  understand  these 
anomalous  phenomena  from  a  model  of  exciton-phonon  coupling  system  with 
inhomogenous  broadening. 

THE  MODEL 

For  inhomogeneously  broadened  materials  such  as  a  PD.-\-LBCML'  film  vhici'i 
contains  linear  chains  of  various  lengths,  there  is  a  distribution  of  the 
exciton  resonance.  The  situation  is  reflected  in  the  observed  absorption 

spectrum  shown  in  Fig.  1(a), ^  where  the  pump  frequency  is  indicated  by  an 
arrow.  It  is  clearly  seen  from  the  curve  that  the  exciton  resonance  at  2.35 
eV  (-527.5  nm)  and  its  vibrational  sideband  at  2.60  eV  are  inhomogeneously 
broadened.  Ue  assume  a  Gaussian  distribution  for  exciton  energies  in  the 
s  amp 1 e , 


g(w)  - 


(1) 


1 

-  exp 

Jlv  tU) 


with  the  parameters  and  Aoj  to  be  determined  by  fitting  the  data.  The 

imaginary  part  of  the  linear  susceptibility  is  calculated^  by  adjusting 
these  parmeters.  As  it  turns  out,  we  can  reproduce  the  experimental 
spectrum  by  choosing  the  most  probable  frequenccy  -  2.35  eV  and  the  width 

Aw  -  0.12  eV.  The  result  is  shown  in  Fig.  1(b). 

With  the  distribution  (1)  in  mind,  we  proceed  to  calculate  the  nonlinear 
optical  susceptibility  due  to  excitions  of  frequency  and  average  our 

result  over  this  distribuiton .  Thus,  we  start  with  the  Hamiltonian 


T  ^  ^  ^ 

H(w  )  —  (j.a'a  +  wb '  b  +  Aa  a(b'+b) 


-  [/i  £(r.p  e 


.i.  iw^t 

+  e  +  h.c. 


(2) 


where  a'  (a)  and  b'  (b)  stand  for  the  creation  (annihilation)  operators  for 
the  exciton  and  the  phonon  with  corresponding  frequencies  and  w, 

respectively.  The  coupling  constant  for  the  exciton-phonon  interaction  is 
denoted  by  A,  and  the  dipole  moment  matrix  element  of  the  exciton  is  given 
by  n.  The  pump  (probe)  field  has  the  amplitude  (E^)  and  frequency 

Wp  (w^)  ,  Only  one  single  phonon  mode  in  the  C^C  stretch  which  couples 

most  strongly  to  the  exciton  is  consiered  here  because  no  qualitative  change 
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is  expected  by  including  more  modes  in  similar  problems. 


Since  we  are  not  interested  in  any  quantity  that  is  sensitive  to  the 
quantum  number  counting,  we  neglect  quantum  fluctuations  and  deal  with  mean 
values  of  the  operators.  Thus,  we  write 

a  -  <a>,  P  -  <b>  (3) 


and  the  equations  of  motion  become  classical.  For  convenience,  we  define 
the  Rabi  frequencies 


n  -  ,uE 

P  .  t  P  ,  c 


and  detunings  , 

A“=cj-a,’  A-c'-w  (5a,b) 

p  P  ^  P  ^ 

0 

In  what  follows,  we  discuss  the  steady-state  solution.^  As  is  always  truc- 
in  pump-probe  experiments,  the  pump  field  is  much  stronger  than  the  probe 

field  which  can  then  be  treated  as  a  small  perturbation.  The  unperturbed 
soluton  can  be  obtained  by  setting  -  0  and  the  results  are^*^ 


Q„  -  n  /(A  -A  n-i7  ) 

0  p'  p  p  X 


(6a) 


/Sq  -  -iAn/(iu;  +  7) 


(6b) 


2  2  2  2 
■where  —  2A  w/(w  +  7  ),  and  the  exciton  number  n  —  ja)  is  determined  by 


2  3  2  2  2  2  2 

A'^n  -  2a  A  n  +  (A  +  7  )  n  -  m  I  -  0 
P  P  P  ^  p  '  p' 


It  is  noted  that  we  have  introduced  phenomenologically  the  damping  rates  7^ 

and  7  for  the  exciton  and  the  phonon  in  the  solutions,  respectively. 

Equation  (7)  may  have  one  or  three  real  roots  depending  on  the 

parameters  A  A  n  and  7  .  There  exists  a  threshold  value  0°  of  the  pump 
p  p  p  X  P 

intensity  below  which  Eq .  (7)  can  have  only  one  real  root.  Our  numerical 

study  shows  that  this  threshold  for  three  real  roots  changes  when  the  other 

parameters  change.  In  fact,  we  find  that  Eq.  (7)  can  have  only  one  real 

root  when  A  <  2  7  no  matter  what  A  and  0  may  be  because  of  the 
P  P  P 

complicated  nonlinear  interactions  between  the  driving  field  and  the 
excitons.  Corrections  5a  and  5/3  due  to  the  presence  of  the  weak  test  field 

are  known  from  numerical  analysis^^  to  be  oscillating  with  frequency  A^ 

around  and  ,  respectively.  Thus  we  have 


iA^t  -iAj.t 

a  —  Q  +  50,6  +  5q.  e 

o  1  2 


iA  t  -iA  t 

^  +  fi/Jj  e  ^  +  5^2^ 


The  susceptibility  experienced  by  the  test  field  is  given  by 
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RESULTS  AND  DISCUSSION 

Equation  (11)  is  employed  to  calculate  the  absorption  spectrum  in  the 
neighborhood  of  -  w  where  the  behavior  of  a  hole  has  been  carefully 


investigated  experimentally.  Parameters  used  in  our  numerical  work  are 

chosen  according  to  experiments.^  They  are  -  2.35  eV,  Aw  -  0.12  eV, 

w  -  0.25  eV.  w  -  2.1  eV,  7  -  O.OOOA  eV,  7  -  0.0002  eV,  A  -  0.01  eV  and 

P  ^ 

n^|/i|^  -  2.5  X  10 '^/x  eV.  In  Figs.  2  and  3,  the  width  and  depth  of  the  hole 

are  plotted  as  functions  of  the  pump  field  intensity,  respectively.  It  is 
observed  from  Fig.  2  that  the  full  width  at  half  the  maximxim  (FUHM)  is 
oscillatory  as  0^  increases.  The  nonlienar  behavior  of  the  hole  depth  is 

also  evidently  seen  in  Fig.  (3).  For  the  range  of  the  incident  laser 
intensity  0.0001  <0  <  0.0003  eV  for  which  measurements  are  made,  both  the 

y  -  p  - 

width  and  the  depth  are  found  in  qualitative  agreement  with  experiments.  We 

have  also  found  that  the  absorption  reaches  saturation  at  about  w  -  0.0011 

P 

eV  not  shown  in  the  figure. 

It  is  perhaps  enlightening  to  point  out  that  the  detuning  is  actually 


a  variable  in  this  model  because  w  is  not  fixed  but  is  characterized  by  the 

X 


distribution  (1).  As  a  consequence,  it  is  always  possible  to  find  for  any 
given  ,  a  Ap  >  2-y^  so  that  Eq .  (7)  has  three  distinct  real  roots.  Since 

different  exciton  numbers  n  imply  different  lattice  relaxation  energies 
which  then  lead  to  different  positions  and  widths  of  the  hole.  Thus,  there 
are  generally  three  distinguishable  susceptibilities  that  have  to  be 
averaged  before  the  integration  over  w^.  In  fact,  our  numerical 


calculations  show  little  or  no  anomaly  when  A  <27  or  n  < 

■'  p  —  X  p 


It  should 


also  be  emphasized  that  extreme  care  should  be  taken  in  the  evaluation  of 


the  integral  in  Eq.  (11).  We  use  10000  mesh  points  for  the  integration 
range  of  2  eV,  about  twice  the  distribution  shown  in  Fig.  1(b) . 

Finally,  we  note  that  the  calculated  hole  width  increases  monotonically 
with  increasing  pump  intensity  for  a  fixed  exciton  frequency  It  is 

therefore  the  inhomogeneous  distribution  (1)  that  causes  the  oscillatory 
variation.  This  implies  that  the  interation  of  the  pump  field  with  excitons 
(both  real  and  virtual)  of  various  frequencies  interferes  and  results  in  the 
width  change  as  well  as  the  saturation  of  the  depth.  To  understand  fully 
the  physical  processes  that  arc  responsible  for  this  anomaly,  a  more  careful 
study  including  dispersion  and  dephasing  effects  are  being  carried  out  and 
results  will  be  reported  elsewhere  in  the  near  future. 
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FIGURE  CAPTIONS 


1(a) .  Optical  density  versus  the  incident  laser  wavelength  observed 
experimentally  for  a  PDA-4RCMU  film  from  Ref.  6. 

(b) .  Absorptive  part  of  the  linear  susceptibility  as  a  function  of  the 
incident  laser  wavelength  calculated  with  parameters  2.35  eV  and 


2. 


Aw  -  0.12  eV. 

FVUM  of  the  hole  in  the  vicinity  of  cj  =  w  in  the  differential 

t  p 

absorption  spectrum  calcualted  as  a  function  of  the  Pv.abi  frequency  of 
the  pump  field. 

Calculated  hole  depth  near  w  ~  w  in  the  differential  absorption 

c  p 


spectrum  as  a  function  of  the  Rabi  frequency  of  the  puir^p  field. 


